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This work is an extended version of the paper [Ij, in which the main 
results were announced. We consider certain classical diffusion process 
for a wave function on the phase space. It is shown that at the time 
of order 10~^^ sec this process converges to a process considered by 
quantum mechanics and described by the Schrodinger equation. This 
model studies the probability distributions in the phase space corre- 
sponding to the wave functions of quantum mechanics. We estimate 
the parameters of the model using the Lamb-Retherford experimental 
data on shift in the spectrum of hydrogen atom and the assumption 
on the heat reason of the considered diffusion process. 

In the paper it is shown that the quantum mechanical description 
of the processes can arise as an approximate description of more exact 
models. For the model considered in this paper, this approximation 
arises when the Hamilton function changes slowly under deviations of 
coordinates, momenta, and time on intervals whose length is of order 
determined by the Planck constant and by the diffusion intensities. 
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1 Introduction 

In this paper we propose a model which, on the one hand, allows one to 
estimate the probability distribution of a quantum particle in the phase space 
in the low temperature heat field. For the first time this problem was solved 
by Wigner |3], but he constructed "quasi-distributions" on the phase space 
which can be negative and hence have no physical sense. On the other 
hand, the proposed model yields one more construction of quantization of 
mechanical systems and can be used in the new approach to foundation 
of the classical quantization procedure. This is an old problem. Various 
approaches to this problem, in particular probabilistic ones, can be found in 
[H El El [TJ [HI [9]. These works essentially infiuenced the author during the 
construction of the present model. 

In this paper we consider the classical model of a diffusion process for 
a wave (complex valued) function on the phase space. The analysis of the 
differential equation of the model shows that the motion in the model splits 
into rapid and slow motions. The result of the rapid motion is that the 
system, starting from an arbitrary wave function on the phase space, goes 
to a function belonging to certain distinguished subspace. The elements of 
this subspace are parameterized by the wave functions depending only on 
the coordinates. The slow motion along the subspace is described by the 
Schro dinger equation. 

Using the assumptions on the heat reason of the diffusions and the corre- 
spondence of the consequences of the model with the known physical experi- 
ments of Lamb-Retherford [10] (the Lamb shift in the spectrum of hydrogen 
atom), we estimate the diffusion coefficients and the time of the transition 
process from the classical description, in which the Heisenberg indeterminacy 
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principle does not hold, to the quantum description in which the Heisen- 
berg principle already holds. The time of the transition process has order 
1/T ■ 10~^^sec, where T is the temperature of the medium. 

The results of this work have been announced in the paper [jj. Proofs of 
theorems 4 and 5 are instructive but rather technical, hence they are exposed 
in the Appendix. The estimate of the parameters of the model is also exposed 
in the Appendix. 

The author is grateful to professor G. L. Litvinov, who was attentive 
to this work, made a lot of editorial comments, and stimulated an essential 
revision of the text, and to professor A. V. Stoyanovsky, who translated this 
paper to English. 

2 Description of the model 

We consider a mathematical model of a process whose state at each mo- 
ment of time is given by a wave function, which is a complex valued function 
(p{x,p), where {x,p) G -R^"", and n is the dimension of configuration space. 
In contrast to quantum mechanics, where the wave function depends only on 
coordinates or only on momenta, in our case the wave function depends both 
on coordinates and on momenta. As in quantum mechanics, it is assumed 
that for the wave functions the superposition principle holds, and the proba- 
bility density p{x,p) on the phase space, corresponding to the wave function 
ip{x,p), is given by the standard formula 

p{x,p) = ip*{x,p)ip{x,p) = (1) 

In the present work we consider the classical model of diffusion process 
for the wave function ip{x,p) on the phase space. It is assumed that each 
complex vector of the wave function is simultaneously in 4 motions: 

the base point of the complex vector moves along the classical trajectory 
given by the Hamilton function H{x,p); 

the base point of the vector moves randomly with respect to coordinates 
and momenta, being in diffusion process with constant diffusion coefficients 

and 6^ with respect to coordinates and momenta, respectively; 

the base point of each vector moves along a random trajectory as a result 
of motions described in the two preceding points, and the vector itself rotates 
with very large constant angular velocity u = mc^/h in the coordinate system 
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related with this point, where m is the mass of the particle, c is the light 
velocity, h is the Planck constant; 

the length of all complex vectors of the wave function at the moment t 
of time is multiplied by exp{abnt/h) (this is a purely technical requirement 
which does not affect the relative probabilities of position of the particle in 
the phase space). 

It is assumed that the wave vector ip[x,p,t) at the point {x,p) at the 
moment t of time equals, by the superposition principle, to the sum of wave 
vectors given by the distribution of vectors (f>{x,p,0) at the initial moment 
of time which get to the point {x,p) at the moment t due to the motions 
described above. 

3 The main results 

3.1 The mathematical model of the process 

Consider the diffusion process on the phase space in which the wave function 
(p{x,p,t) at the moment t satisfies the differential equation 

where A,,^ = gl^— -^)^ + h^ E + — ^, 0) 

where H{x,p) is the Hamilton function; and are the diffusion coefficients 
with respect to coordinates and momenta, respectively. 

If we omit the last summand in equation ([2]), then we obtain a first order 
partial differential equation dip/dt = A(f, where 

A _y(9H dip OH dip \ t / " dH\ 
ir[\dxkdpk dpkdxk) h\ ^^dpk^^)'^' 

This part of equation ([2]) describes the deterministic component of the 
motion of complex vectors ip{x,p, t) along the characteristics of the equation. 
According to the equation, in this motion the base point of each vector 
moves along the classical trajectory given by the Hamiltonian H[x,p), and 
the vector itself rotates at each point of the trajectory with the angular 

velocity u;' = li^H-EU§-^Pk). 
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Note that in the case when the configuration space is three dimensional 

and H = E = c\Jm^c^ + p^, we have 

where r = mc^dt/H, in accordance with the formulas of special relativity 
theory, is the proper time in the coordinate system related with the particle 
moving with the momentum p. I. e. in this case, the vector whose base 
point moves along the classical trajectory, rotates with the constant angular 
velocity u = mc^/h in the coordinate system related with this point. 

On the contrary, if in the right hand side of equation ([2]) we leave only 
the last summand of the form ([3]), then we obtain the equation 

dip 2^f 9 ipkV ,,2^*92 abn 

This equation describes the diffusion component of the motion of vectors 
ip{x,p,t) in the phase space. In this motion, the base points of the vectors 
move according to the classical homogeneous diffusion process with the dif- 
fusion coefficients with respect to coordinates and momenta equal to and 
6^, respectively. And the vector itself is parallel transported during small 
random transports from a point {x,p) to the point {x + dx,p + dp), and 
its length at moment t is multiplied by exp{abnt/h). Note that the parallel 
transport of vectors on the phase space is given by a connection expressed 
by the following formula: L(^dx,dp)V{^^P) ~vi.^i'P) ~ —{i/fi)ip{x,p)pdx, where 
L{dx,dp)'p{x,p) is the parallel transport of the vector (p{x,p) from the point 
{x,p) along the infinitely small vector {dx,dp). 

In the particular case when the configuration space is three dimensional, 
such the connection on the phase space is related to synchronization of mov- 
ing clocks at the points of the phase space. Indeed, if a particle with coordi- 
nates X = {xi, X2; X3) moves with the velocity v = {vi, V2] v^), then, according 
to formulas of special relativity theory, proper time is expressed through the 
observer time t by the formula 

t — xv/c^ 
^1 - v^/c^ 

where xv = XiVi + X2V2 + X3V3 is the scalar product of the vectors x and v, 
and c is the velocity of light. 
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(6) 



For a free particle with the momentum p = {pi] P2'-, Pz) and the stationary 
mass m, the energy E = c\/p^ + vri^c? and, respectively, 



V = /l - — = (7) 

V^p2 _|_ 77^2 c2 ' y -|- TTJ^C^ 

Substituting these expressions into ([6]), after computations we obtain: 

T = ^1^. (8) 

Consider the distribution of complex vectors (p{x,p, t) on the phase space, 
rotating with constant angular velocity u = mc^/h in the proper time r, 
which, at the moment r = 0, are equal to one and the same vector (pQ. We 
have 

o 

V?(x, p, t) = ipo exp f -z-^r j . (9) 

Substituting formula (IE]) into this expression, we obtain 

/ \ / —iiEt — xp)\ 

= v^oexpl 1. (10) 

Hence, if 0)^^(2;, P, ^) is the shift of the vector y9(x,p, t) = {poexp{—imc^T/h) 
by the vector Ax along coordinates x without change of proper time, then 

L{Axfi)'f[x,p,t) = Lp{x,p,t)expl — - — 1. (11) 

In the limit of infinitely small Ax we obtain the required formula for this 
case: L(^dxfi)'^{x,p) — ip{x,p) = —{i/h)ip{x,p)pdx, where dx is the infinitely 
small shift of coordinates x. 

On the other hand, if we have a shift -^v(o,Ap) of the vector {p{x,p,t) = 
ipQexpl—imc^r/h) along momentum by Ap without change of proper time, 
then, since in the special relativity approximation, acceleration does not 
change the proper time of a particle, we have the equality L(o,Ap)V^(a^?P5 ^) = 
ip{x,p,t) and L(Q^dp)^{x,p) - ip{x,p) = 0. 

Hence, by linearity of L(^dx,dp) with respect to {dx,dp), we obtain the re- 
quired equality in the general case: L(^dx,dp)'^{x,p)—(p{x,p) = —{i/h)ip{x,p)pdx. 

Note that under these assumptions, the derivation with respect to the 
vector of infinitely small shift along the k-th coordinate corresponds to the 
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differential operator Dx^ = d/dxk — ipk/^, and the derivation with respect 
to the shift along the k-th momentum corresponds to the usual differential 
operator Dpj^ = d/dpk, where k = 1, ...,n. 

Note also that these operators of shift along coordinates and momenta 
do not commute. The commutators of these differential operators read 



—i/h and 



0, where k ^ j and k,j = 1, 



n. 



Thus, the shifts along coordinates and momenta of wave functions on the 
phase space realize a representation of the Heisenberg group. 



3.2 Analysis of the diffusion component of the equa- 
tion 

Consider the diffusion equation ([5]) in more detail. This equation can be 
represented as follows: 

dip 2 a 7 2a ^'^^ 

— = Aa,b'P = a A^v9 + b ApV9 + — V^, (12) 

where 

n n o2 

k=l k=l ^Pk 

It is natural to call the operator Aa^ by the diffusion operator for the 
representation of the Heisenberg group with the diffusion intensities a and b 
with respect to coordinates and momenta, respectively. 

Let us look for a solution of equation (fT2i) in the form 

(fi{x, p, t) = ip^{x, p, t) exp {^-^^ ■ (13) 

Substituting this expression into equation f|T2l) . dividing both parts of the 
equality by exp{ixp/h) and transferring the summand {nab/h)ip^ into the left 
hand side, we obtain the equation 



a/ nab " /2'9V° ,2/^V° 2ixfc(9/ 
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where = ip^{x,p^t) is some function. 

To solve equation (fT^ . let us decompose the function ip^{x^p,t) into the 
Fourier integral with respect to p, i. e., let us represent the function t) 
in the form 

1 /" / 2 yp \ 

where = J^^^(fi%x,p,t) = ^^^^^^^^ j^^(p^{x,p,t) expi^—jdp. 

(16) 

Substituting this expression for ip^{x,p,t) into the equation (JT4l) . we obtain 
that ip^{x,y,t) satisfies the equation 

-V = i:fa^?^-!^^^^^V'») (17) 



dt h dxl h 

The right hand side of this equation is a self adjoint operator with discrete 
spectrum consisting of negative numbers. 

Indeed, the equation for eigenvalues of this operator reads 



where x = xi^iU) is a function of x and y. This equation is the stationary 
Schrodinger equation for harmonic oscillator, and it is well studied (see, for 
instance, [13], p. 94). 

In particular, it is known that on the set of functions which tend to zero as 
X tends to infinity, equation (fT8|) has discrete spectrum consisting of negative 
eigenvalues Ai > A2 > . . .. The greatest eigenvalue Ai = —nah/Ti corresponds 
to the eigenfunction Xii^iU) = (b/anh)^^^ exp (—b(x — y^ /2ah). The next 
eigenvalues are less than Ai, and the difference is greater than or equal to 
ab/h. 

Since the eigenfunctions Xk{x,y) of the operator (|T8l) form a complete 
system of functions in the class of functions tending to zero as x tends to 
infinity, an arbitrary function tp^{x,y,t) from this class can be represented 
as a series 

00 

V'°(x,?/,t) = J2(^k{y,t)xk{x,y), 

k=l 
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where Ck{y,t) = tlj^{x,y,t)xk{x,y)dx (19) 

are the coefficients of the decomposition of the function tlj^{x,y,t) with re- 
spect to eigenf unctions Xk{x,y). 

Substituting the expression of the function ip'^{x,y,t) in the form of this 
series into equation (fTTl) . we obtain that this equation in the orthonormal 
basis of eigenfunctions Xk{x,y), k = 1,2, .. ., sphts into an infinite system of 
equations: 

dck{y,t) ( nah\ 710 
— di — = + — \ck[y,t) k = l,2,... 

where Ai + nab/h = and + nab/h < —ab/h for k > 1. 

Hence ci{y,t) = ci{y,0), and Ck{y,t) = Ck{y,0) exp{{Xk + nab/h)t) expo- 
nentially decay with time for k > 1. Hence the summand in corresponding 
to the ffist eigenvalue will give the main contribution into the function ip'^ 
after time of order {fi/ab). 

Thus, we have obtained that a solution of equation (ITTl) . after time t of 
order h/ab, becomes exponentially close to the function ip'^{x,y), where 

'^^{.x,y) = lim i)^{x,y,t) = ci{y,0)xiix,y). (20) 

t — ^00 

Respectively, since by definition ip^{x,p, t) = J^ntp^^x, y, t) and Fourier trans- 
form is continuous, we have 

^\x,p) :Fni^\x,y) = \im :Fni^\x,y,t) = rn{c^{y,Q)Xi{x,y)). (21) 

t — >OQ 

Since we will not use other eigenfunctions, introduce the notation 

X{x - y) =^ xiix,y) = (b/anhr^^exp (-6(x - yy/2ah) . (22) 

To make the formulas shorter, let us also denote il!{y) '= ci(?/,0), where 
ciiy, 0) is expressed by the formula (|T9|) with k = 1 and t = 0. That is, 

i^{y) = {x, y, 0)xi{x, y)dx = ij^{x,y,0)x{x - y)dx. (23) 

Jr" Jr" 

Since by formula (IT^ ip{x,p,t) = ip^{x,p,t) exp{ixp/h), then by formulas 

f l2T|) . fl22l) and using the notation tp{y) Ci{y,0) above and also the equal- 
ity (1251) and notation flTB]) . we obtain the following theorem. 
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Theorem 1. Let !f{x,p,0) be an arbitrary function such that Fourier 
transform of the function (f{x,p,0)exp{—ixp/h) with respect to p tends to 
zero as X — > oo. Then the solution ip{x,p,t) of the diffusion equation ^ 
exponentially with time (with the number in the exponent equal to —abt/h) 
tends to a stationary solution of the form 

M^,P) = ^lim = Jmxix - y)e-'^y-^^^'Hy, (24) 



where ip{y) = j^—^ j ip{x,p,0)e''^y ""^P^^ix ~ y)dpdx, (25) 



and x{x-y) = {^] e-''("-^)'/(2aft)_ ^26) 

\aTrn J 

Note that x^{^ ~ u) is the probability density of the normal distribu- 
tion with respect to x with the mathematical expectation y and dispersion 
ah /{2b). If the quantity ah /{2b) is small, then the function x^{.^~y) is close 
to the delta function oi x — y. 

The composition of expressions ( l25l) and ( IMl) yields a projector Pq form 
the space of wave functions defined on the phase space onto certain subspace. 
The elements of this subspace are parameterized by functions of the form 
'?/'(?/), where y G -R", i. e., by wave functions on the configuration space. 

Theorem 2. The projection operator Pq given by composition of expres- 
sions (E^ and [24\ ), has the form 

Pof = / f{x,P)e 4afi e e 2^ dx dp . (27) 

The operator Pq is self-adjoint and commutes with the operator Aa^b- 

Proof of this theorem is obtained by substitution of formulas ( l25l) and 
f l26|) into fl2^ . by an algebraic transformation of the number in the expo- 
nent, and by computation of the integral over y. The integral over y is the 
Fourier transform of the exponent of a quadratic polynomial, whose analyt- 
ical expression is known. The author performed the computations using the 
system Mathematica [H], supporting symbolic mathematical computations. 
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The commutativity of the operators Pq and A^,;, follows from the fact that 
the orthogonal projector Pq distinguishes the subspace of eigenvectors of the 
self-adjoint operator Aa^t with the zero eigenvalue. 

Formulas fl24|) and ([1]) imply the following statement. 

Theorem 3. If iplx) is a wave function on the configuration space and 
ip{x,p) is the wave function on the phase space corresponding to it by for- 
mula ( P^ , then the probability density in the phase space p{x,p) = 
is expressed by the formula 

In contrast to quasi-distributions 

1 r ( i*f f ix'p^ 



dx' 



defined by Wigner [3], the density p{x,p) in the phase space, given by expres- 
sion fl28l) . is always nonnegative. Its expression differs from the expression 
for the Wigner function by exponents under integral, which give the densities 
of distributions close to delta functions. 

To prove Theorem 3, one should substitute into formula ([T]) expres- 
sion (1211). We obtain 



pix,p) = = — i-^ I ^{y)riy')x{x - y)x{x - y')e-'^y-y'^P/^dydy', 

(29) 

where x{.x — y) is given by relation (1261) . After substitution of (126|) into (!29|l . 
the change of coordinates y = x + {x" — x') /2 and y' = x + {x" + x')/2 under 
integral and a transformation of the numbers in the exponents, we obtain 
formula (|28l). 

The algebra of observables given by real functions on the phase space but 
averaged over the probability densities of the form (l28l) . has been studied 
in [ig. 

The function p{x) of density of probability distribution in the configura- 
tion space is expressed through the density p{x,p) in the phase space by the 
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formula p(x) = Jj^n p{x,p)dp. Hence, integrating expression (12^ over p, we 
obtain the following statement. 

Corollary 1. Ifiplx) is a wave function on the configuration space, then 
the corresponding probability density p{x) in the configuration space is given 
by the formula 



where xi^ ~ u) is given by [2D\) . That is, p{x) is obtained from \iIj{x)\'^ 
by smoothing (convolution) with the density of normal distribution with dis- 
persion ah /{2b), and the exactness of defining coordinate is bounded by the 



As it is known |16], in quantum electrodynamics the minimal error of 
measuring coordinates of an electron in the stationary system is bounded by 
the quantity h/{mc), where m is the mass of the electron, and c is the light 
velocity. Hence the statement of Corollary 1, although not corresponding to 
non-relativistic quantum mechanics (in which it is assumed that coordinates 
can be measured with any degree of exactness), does not contradict with a 
more exact theory, quantum electrodynamics. 

If one assumes that the diffusion is induced by heat action on the elec- 
tron, then the diffusion coefficients with respect to coordinates and momenta 
are expressed in statistical physics (see, for example, [T7], Ch.7, §4 and §9) 
through the temperature T by the formulas = kT/{m'~f) and 6^ = 
^kTm, where k is the Boltzmann constant, m is the mass of the elec- 
tron, 7 is the friction coefficient of the medium per unit of mass. Hence 
a/b = (7m)~^ and ab = kT. That is, in this case, the quantity a/b, which 
enters expression fl2Bl) and determines the dispersion of smoothing in Corol- 
lary 1, does not depend on the temperature. On the other hand, the time 
t of the transformation process determined in Theorem 1, has the form 
t ~ h/{ab) = h/{kT) = • 7.638 ■ IQ-^^ ^ ^ore detailed formulas for 
the estimate of the quantity h/{ab) are given in Appendix 3. 

3.3 Analysis of the model of the process 

Let us return to the study of the main equation ([2]). Taking into account 
the estimate made at the end of the previous subsection, let us consider the 
quantity h/ {ab) in equation ([2]) as a small parameter, and let us assume that 
coordinates and momenta change a little at this time in the classical motion 




(30) 



quantity 
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defined by the Hamiltonian H{x,p), and also let us assume that the function 
H{x,p) and all its derivatives grow at infinity no faster than a polynomial. 

Theorem 4. The motion described by equation ^ asymptotically splits 
as h/{ab) into a rapid motion and a slow one. 

1) As a result of rapid motion, an arbitrary wave function 99(0;, p,0) turns 
at the time of order h/{ab) into a function of the form [24\): 

V{x,p) = j^J^ [mx{x - y)e-'-<y-^y^/Hy, (31) 



n/4 



where x{x~y)=\-^\ e-''("-^)'/(2aa)_ ^33) 

The wave functions of the form ( fjlj) form a linear subspace. Elements of 
this subspace are parameterized by wave functions ip{y) depending only on 
coordinates y G -R". 

2) The slow motion starting from a nonzero wave function (f{x,p, 0) of the 
form l[3l\) from this subspace, goes inside this subspace, and is parameterized 
by the wave function ip{y, t) depending on time. The function iply, t) satisfies 
the Schrodinger equation of the form ihdip/dt = Hip, where 

^X{.x - y)xix - y')err^y-y''^Pip{y\ t)dy'dxdp, 



and x(x — y) is given by formula [3^) . 

Proof of Part 1 of Theorem 4 is postponed till Appendix, due to its large 
volume and technicalities. 

Proof of Part 2 of Theorem 4. In the first Part of Theorem 4 it is 
stated that after the rapid motion, the initial distribution ip{x,p, 0) turns to 
the form with just a small difference with ( l3Ti) . After the slow motion the 
distribution remains in the class of functions of the form (13T!) . but changes 
in time. 

To study the slow motion, let us look for a solution of equation ([2]) in the 
form ( |3T1) in which ip = ilj{y',t) is considered as dependent on time. 

Let us substitute expression (13T1) into equation ([2]). Since by construction, 
this expression for (p satisfies equation Aafi(p = 0, then after substitution of 
expression fl3Tl) into equation ([2]) and after dividing both parts of the equation 
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by exp (ixp/h), we obtain the following equation: 
1 f d^{y',t) ( iy'p\ , 



(27r;i)"/2 



k=l 



If in the obtained equation one opens the brackets, reduces the similar terms, 
multiplies both parts of the equation by l/(27r^)"'/^x(x — y) exp(iyp/h) and 
integrates by p and by x, then, taking into account the equality J^„ x^{^ ~ 
y)dx = 1, we obtain the following equation: 

- = --i7V^ or ^n^ = H^, (35) 



where the operator H is obtained from the function H by formula (!33|) . 
required in Part 2 of Theorem 4. 

Theorem 5. If ^ is a small quantity and H{x,p) = ^ + V{x), then 
the operator H , up to terms of order ah/b, has the form 

Proof of Theorem 5 is given in Appendix 2. 

The first two summands in formula (!36l) give the standard Hamilton op- 
erator. The last summand is a constant and can be neglected. The previous 
summand before the last one will be considered (due to the smallness of ah/h) 
perturbation of the Hamilton operator. 

Assuming that the deviations in the spectrum of the hydrogen atom (the 
Lamb shift) observed in the Lamb-Retherford experiments [10] are induced 
by the previous to the last summand in formula fl36|) . one can estimate the 
quantity a/h. The computations by the standard method of perturbation 
theory analogous to the computations of [TH], give the following estimate: 
a/h = 3.41 ■ lO'^sec/g (see the computations in Appendix 3). Hence, the stan- 
dard deviation of the normal distribution x^, with which we make smoothing 



in formula fl30|) . has the form ^ah/{2b) = 4.24 ■ lO^^^cm. This quantity 
is much less than the radius of the hydrogen atom, and it is close to the 
Compton wave length of the electron h/{mc) = 3.86 ■ 10~^^cm. 
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4 Conclusion 



In this paper it is shown that the standard quantum mechanical description 
of a process can arise as an approximation of certain classical model for a 
diffusion process for a wave function in the phase space. The computations 
show that the proposed model in the form of differential equation ([2]) de- 
scribes physical processes in a rather adequate manner in the standard cases 
for standard Hamiltonian. But this model can be applied as well for com- 
putations of processes with a nonstandard Hamiltonian or a Hamiltonian 
rapidly changing in time, as in sudden perturbations [19j or for periodically 
changing potential with frequency of order ah/h, and it can be compared 
with experimental data. 
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APPENDICES 



Appendix 1 



Proof of Part 1 of Theorem 4 

Let ip{x,p,t) be a solution of equation ([2]). In the notations ([3]) and (jlj) 
equation ([2]) takes the form 

^ = Aip + Aa,b^, (37) 

where A is a skew Hermitian operator, and Aa^b is the self-adjoint diffusion 
operator. 

Consider the derivative with respect to time of Hv^lP (the square of the 
norm of the function <f). We have 

^^'^''^^ + = {Av + Km^) + {^^,^^ + ^aM■ 



dt dt 'df^' dt 

Hence, by linearity of the scalar product ( ; ) and by the equalities {Aip; ip) = 
— {(p; Aip) and {Aa^^V^ Lp) = {^p; Aa^fp), we obtain the equality 

2|MI^ = ^^ = 2(A,,^;^) or ^M.^(A„,»,;»,). (38) 

def 

Denote by = the normalized function (p. If the function (p 

satisfies equation (1371) . then, taking into account formula (1381) . we have 

dip> d f (p \ 1 dip ip d\\(p\\ 



dt dt \\MJ \M dt ||(^||2 dt 



Aip+Aa^b^-{Aa,b^;(p)ip. (39) 



By definition of the function (p its norm ||<^|| = 1, and it is natural to call 
equation ( |39l) the equation for the normalized wave function. 

Consider now the projection of the function ip{x,p,t) on the subspace of 
stationary solutions of the diffusion equation ([5]). 

Let us represent the function ip{x, p, t) in the form ip = yjo+V'i, where ipo = 
Poip, ipi = {l—Po)'p = Pi^p, and Pq is the projector from Theorem 2. The self- 
adjointness of the projection operator Pq implies the equalities {(po; ipi) = 
and llv^olP + llv'ilP = 
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For the proof of the first Part of Theorem 4 one needs to show that the 
quantity Hv^ilP = 1 — Hv^olP becomes small after time of order h/{ab). 
To make formulas shorter, introduce the notation 

n{t)=^ \\Mt)\\^ where = Po^. (40) 



Statement 1. If(p satisfies equation ^^3l\) , then the quantityqit) = Hv^oH 
satisfies the differential equation 

f, = a{t)^l^^ + - r/)r7, (41) 

where (3{t) > Pmin = 2ab/h and \a(t)\ < amax = 2 ■ max^^ \ \Aifo - -Po-4<^o||- 
The maximum in the latter expression is taken over all normalized functions 
(fo from the subspace of stationary solutions of the diffusion equation l^, for 
which the function {(fo^Xjp)]"^ gives a probability distribution in the physical 
region of the phase space for the given process. 

To prove Statement 1, let us consider how the quantity Hv^olP changes in 
time, where ipo = Pqi^ and (p satisfies equation (1371) . 

We have the following equalities which follow from the definition of the 
function ipo, from independence of the operator Pq of time, from equality (|39|) . 
from linearity of the scalar product ( ; ) with respect to each argument, from 
commutativity of Pq with Aa,?,, from self-adjointness of the operators Pq and 
Aafi, from skew Hermitian property of the operator A and from the projection 
property Pq = Pq'. 

(substitute expression for = ^jq + V^i, use 
the linearity of operators and equalities A* = —A, Aa,b^o = 0) 
= -(v^i; ^Vo) - {Aipo; Vi) - 2(Aa,feV5i; v^i) ■ (v^o; V^o) 

/. V'o \ 




IV^oll IIV^il 



■ ll<^o| 



llv'ill W^iW/ 

-{{<fi;A<fo) + {A(po; <^i))||v?i|| \\(po\\-'2{Aa,b<fi; <^i) | Iv^il H Iv^ol (42) 
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where ipi and <^o are the normahzed functions y^i and ip^, respectively. 
Since rj = \ \ipo\ p and \ \(po\ P + I Iv^il P = I P = then | |v?i| | = y/T^-rj. 
Substituting these equahties into equation fH2l) and introducing notations 

def def 

"(^) = - ((<^i;^<^o) + (^<^o;<^i)) and/3(t) = -2(A„,b(^i; (^i), 

we obtain the equation from Statement 1: rj = a{t)y/l — ri^yrj + P{t){l — rf)rj. 

For the estimate of note that since ((^i; Aipo) = {Ai^q; <^i)*, where 

at the end of the equahty the star * denotes complex conjugation, then 
a{t) = 2 ■ Re{{—Ai^Q]i^i)), where Re denotes the real part of a complex 
number. Hence we obtain the following estimates: 

\a{t)\=2\Re{{-A^o;^i))\ < 2\{A^o; = 2\{A^o - PqA^o] = 
2|((A-Po^)v^o;v^i)| <2||(A-Po^)<^o|| ||(^i|| = 2||A<^o-i^o^<^o||. (43) 

The latter inequality in the previous formula is obtained using the Cauchy- 
Schwartz inequality estimating the absolute value of the scalar product. The 
function —PqAi^q, which is orthogonal to the function ipi, has been added to 
Aipo, in order to obtain the projection of the function AipQ onto the subspace 
orthogonal to the space of stationary solutions of the diffusion equation ([5]). 

Let us now estimate the quantity P(t)= — 2(Aa^fe(^i; (^i). To this end, 
decompose the function (fi with respect to the orthonormal eigenfunctions 
Xi of the self-adjoint operator Aa,?,. We have (fi = J^iZo^^iXt- Since the 
vector (fi is orthogonal to the vectors with eigenvalue equal to 0, we have 
Aa,fe<^i = Ei: x^^o^iCiXi and Ei: A,^oC? = 1. Hence 

oo 

{Aa,b<fU<fl) = ( J2 ^iCiXhJ^'^iXi) = ^i^i ^ ^max<?i = >^max, (44) 

i: Ai^O i=l i: A^^O i: A^^O 

where Xmax is the maximal nonzero eigenvalue of the self-adjoint operator 

Using this inequality and also knowing that Xmax = —ah/% (see Proof of 
Theorem 2), we obtain the estimate required in Statement 1: 

Statement 1 is proved. 

Statement 2. As h —>■ the estimate of the quantity O-max has the form 
ctmax = (C + Ofi{l))/\/h, where Oh{^) is an infinitely small quantity with 
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respect to h, and the constant C is determined by the maximal value of the 
absolute values of the first and second derivatives of the Hamilton function 
H{x,p), when the coordinates and momenta belong to the physical region of 
values for the given process. 

Proof of Statement 2 will be postponed until the end of this section. 

For the estimate of the time of the transformation process described in 
Theorem 4, let us rewrite the equation from Statement 1 in the following 
form: 

^ <^ ^ 

a{t)./T^^ + P{t){l-r])r] 



a{t)^/l - r]^ + f3{t){l -r])ri J -amaxV'^ - Vy/V + Pminil - V)V 

(45) 

where the integral is taken over the interval (e, 1 — e), on which the de- 
nominator of the expression under the latter integral is positive. Simple 
computations show that for this the following inequalities should hold: 

(46) 

If one takes into account that (3min = 2ab/h (see Statement 1) and amax = 
{C + on{l))/Vh (see Statement 2), then, ash ^ 0, the right hand side of the 
previous inequality can be represented in the following form by decomposing 
into the Taylor series: 




/ '^max _ ^max Q f '^"1"^: \ ^ _ C H _|_ ^ / ^ 



2 y 4 (^min Pmin \ /3min / Ao^b"^ \oh 

This and inequality fH^ imply that the number e can be chosen to be arbi- 
trary (small, as ^ ^ 0), satisfying the inequality 

1 C^h fh 



The last integral in inequality fl45l) of the form 

^^ %f Y , (48) 
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has been computed using the system Mathematica 5.0 (see [S]), and the 
result was decomposed into the series with respect to Prnin as Pmin oo. 
These computations yield the following equality: 

^ _ Aarctanhjl -2e) ^ 4(1 - 2e) amax , q ^ \ ^49) 

where arctanh{x) = {ln{l + x) — ln{l — x))/2 is the hyperbolic arctangens. 

If one assumes that e and h are small quantities, then, substituting into 
this equality the expressions for P^in = 2ab/h and amax = (C + ori{l))/Vh 
from Statements 1 and 2, and decomposing the obtained expression into a 
series with respect to e, one can obtain the following estimate: 

t. = (-ln. + o(.))- + ^-AiZ_y +o(^). (50) 

Inequality (145!) . formulas (1471) . (148|) . (l50l) . and the definition of the func- 
tion ri{t) (140|) immediately imply the following statement. 

Statement 3. Let (p{x,p,t) be a normalized solution of equation 

def 

and, at the initial moment for t = 0, let the following inequality hold: ri{0) = 
\\Pq(P{x,p, 0)11"^ > e, where Pq is the projection operator from Theorem 2 and 
e is an arbitrary number satisfying the inequalities 



1 c^h /h 

2 i^^^'l^ 



Then, for any t > ts, where 



h_ jl + oje)) c fhy 



+ + +0^^j, (51) 

the quantity rj{t) =^ ||Po'^(^)P)^)|P > 1 — i- e., the square of the distance 
from the function ip{x,p,t) to the subspace of stationary solutions of the 
diffusion equation, described by Theorem 1, will be less than or equal to e. 

Part 1) of Theorem 4, being proved at this section, immediately follows 
from Statement 3. For conclusion of the proof it remains to prove State- 
ment 2. 

Proof of Statement 2. In the proof of Statement 2 we shall need 
computations of some integrals containing the function xiv) given by ex- 
pression fl26|) from Theorem 1. The results of computations of these integrals 
are listed in the following two Lemmas. 



20 



Lemma 1. Let x be the Fourier transform of the function X; where 
x{y) — (b/anh)^^'^ exp{—by'^ /{2ah)). Then the following equality holds: 

On the contrary, the Fourier transform of x yields x- I- e., 

One has more general integrals: 

(27r^)n/2 X„ - y)e-^dy = xip - k)e-^ , (54) 

X f ik(x~y) ip{x — y) 

{2T:hY/2 - = X{x - y)e-^. (55) 

Functions x o-nd x satisfy the following relations: 

X(a)x(/3) = x(^^^^ anc? x(a)x(/3) = X^^^^ X^^^j- (56) 

The derivatives of the functions xiy) o,nd x{k) have the form 

^ = -b/{ah) y,x{y); ^ = -a/m k.xik)- (57) 

The functions x^ o.iT'd X^ are the densities of the normal distribution in the 
configuration space and in the space of momenta, respectively, with the zero 
mathematical expectations and the dispersions equal to ah /{2b) and bh/{2a). 
I. e., the following equalities for the function x hold: 

I x^{v)dr}^l] / r]iX^{rj)dri^O, i^l,...,n; 

/ ViVjX^i.v)dv = 0, for i j; / ViX^iv)dv = ah/ {2b), i = l,...,n. (58) 

The other moments have order o{h) . Analogous equalities hold for the func- 
tion X- 

ii?" JR" 

I UjX'{Odi = 0, fortj^j; f ^^f{0d^ = bh/{2a),z = l,...,n. (59) 

Jii" JR^ 
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Besides that, the following equalities hold: 

^2-j^f^)n/2 /fj„ Vi^jX{v)x{Oe~drid^ = 0, for i ^ j; (60) 
/ Vj^jXiv)xiO(^^dr]d^ = ^, j = 1, . . . , n. (61) 



(27r;i)"/2 }r 

All integrals from Lemma 1, except for the latter one fl^T]) . are well known. 
The latter integral has been computed by the author using the system Math- 
ematica 5.0 [14j. 

Lemma 2. Let 

D = ^^Jy^r^i, Xiv) x{v7V2) x{aV2) ^{^i)e^i^^^r,i'HiHa2)in_ (g2) 
Then the following equalities hold: 

1) / Ddr] di dri'di' =1; 2) / mDdr] di drj'di' = 0, j = l,...,n 

3) / {.Vj + v'i)Ddr] dC, di]'d^' = 0, j = 1, . . . ,n; 

4) / ^.Ddr] d^ drj'd^' = 0, j = l,...,n; 

5) / {^,+ej)Dd7]d^dr]'de = 0, j = l,...,n; 

6) / miVj + '7!)^c?^ c?^ c?'7'c?^' = 0, z, j = 1, . . . , n; 

7) / (^i + OijDdr] d^ dri'd^' = 0, i, j = 1, . . . , n; 

8) / r]ir]jDdr] d^ drj'd^' = 0, for i ^ j; 

I ri'^Ddri d^ dri'd^' = ah/(2b), j = l,2...,n; 

9) / iv^ + v'i)iVJ+v'J)Ddr] d^ drj'd^' = 0, fori^j; 

/ (rii + ri'.fDdri d^ dri'd^' = ah/(2b), j = l,2...,n; 

10) / ^,^,Ddr]d^dr]'d^' = 0, for i ^ j- 
I (^-Ddr] di didi' = hh/{2a), j = l,2...,n; 



iRi" 
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11) / iC. + m^,+^j)Ddrjd^dr]'de = 0, fort^j- 
J n't" ■' 

f (0 + Q'Ddr] d^ dr^'de = bh/i2a), j = 1, 2 . . . , r 

12) / rfi^jUdr] d^ drj'd^' = 0, z, j = 1, . . . , n; 

13) / ri'^i'pdTidididi' = 2,j = l,...,n; 

14) / r^,{i,+ir)Ddr^didr^'d^' = for t ^ j; 

I vA^j+^',)Ddvd^dr]'dC = th/2, j = l,2...,n; 

/ UVj + v'j)Ddr] d^ d7]'d^' = th/2, j = l,2...,n; 

16) j {S,^ + 0{ri,+ri'^)Dd7]did7]'di' = for i ^ j- 



n. 

Proof of Lemma 2 is based on the use of formulas from Lemma 1. 
Let us compute here the first integral of Lemma 2. The other integrals are 
computed in a similar way. 

For the first integral, after substitution into it, instead of of expres- 
sion (j62l) . we obtain 

= {integrating over ^ by formula fl53|) } 

= (2.;.)'/^2»/^i /(^)^(^^/^)^(^^/^) x(^+V)e^(''-^^'/^^^'/^^r,ciV^e' 
= {integrating over by formula ll53l)} 

= X(^)X (^7^2) X (V2^ + V'/V2) xiv + vVvdv' 
= {apply formula ( l56l) to the two middle factors x} 
x'^{v)x^{v + rj')drjdrj' = 1 {by formula fl58l) }. 



In the proof of Statement 2 we shall also need the averaging of functions 
F{x,p) on the phase space (x,p) G R^"' with respect to the density of the 
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distribution 

K = 7^[ ^{x)r{y)e'^dy. (63) 



This density looks similar to Wigner's quasidistribution, but does not coin- 
cide with it. 

Denote by Fiy^ the average of the function F{x,p) with respect to the 
density W^. That is, 

" (27r\)" Irb^ Fi^^P)^i^)^*iy)(^^'^^dydxdp. (64) 

Lemma 3. If F{x,p) is a smooth function which, together with all its 
derivatives, grows at infinity no faster than a polynomial, and ip{x) is an 
arbitrary smooth complex valued function rapidly decreasing at infinity, then 
the following equality holds: 

lim Fvi// = / F(x,0)^p(x)'ijj*(x)dx. 



Proof of Lemma 3. Let us make the change of variables k = p/h under 
the integral (1641) . and let us integrate over y. We obtain 



= / F{x,hk)^{x)r{y)e'^''-''^'dydxdk 



1 



F{x,hk)ip{x)^lj*{k)e-'''^dxdk, (65) 
where i>*{k) = \ [ i)*{y)e^y^dy (66) 



^— / 

(27r)"/2 jRr 

is the Fourier transform of the function ip*{y). Since the function ip*{y) 
is rapidly decreasing, its Fourier transform ip*{k) is also a function rapidly 
decreasing at infinity. 

Since his a. small quantity, let us decompose the smooth function F(a;, Uk) 
over hk by the Taylor formula up to the terms of first order. We have 

" dF 

F{x, hk) = F{x, 0) + hY, ki — {x, 9hk), 

i=l ^Pi 

where 9 = {9i, ft„) and < 6*4 < 1 for i = 1, ...,n. 
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Let us substitute this expression of the function F{x,hk) into (^^. We 
obtain 



(27r)"/2 jR2n 
1 



h 



F{x,0)^{x)^/j*{k)e-"^dxdk 
" dF 

kij-{x, ehk)^{xW{k)e-'''^dxdk 



(67) 



(27r)"/2 7^2" ^ 

Let us estimate the coefficient before h in the second summand of the 
obtained equahty. We have 



(2vr) 



j T.''^Q-i^^^^f')^i^)'P*W('~''''d^dk 



< 



(27r)"/2 



" dF 



7 — ^ — 777 hm 

(27r)"/2 ''^oo 



Dr={\\{x,k)\\<r} 



< 



< 



1 



(27r)'^/2 
1 



hm 



max 

J (x,k)£Dr 

Dr 



. ,„ hm / max 

(27r)"/2 r^ooj (x,k)£Dr 



\i){x)\ \ip*{k)\dxdk 

\iIj{x)\ \^*{k)\dxdk 
I \%l){x)\ \i)*{k)\dxdk 
\tlj{x)\ \-^*{k)\dxdk 



" dF 

y2ki—{x,ehk) 

" dF 

y2kiT:—{x,9hk) 



" dF 



= \im J 0{x'\k''')\ilj{x)\ \i)*{k)\dxdk = M. 

Dr 

The latter equalities follow from the fact that by the statement of Lemma 3, 
the expression in the integral, standing under the operation max, grows no 
faster than a polynomial of certain degree with respect to each variable, 
and |'?/'(a;)| and \ip*{k) \ are rapidly decreasing functions (decreasing at infinity 
faster than any power), and the limit as r — >■ oo of the integral of a positive 
rapidly decreasing function exists and is equal to some M. 

Formula (1^ and the boundedness of the coefficient before /i — > in the 
second summand of this formula imply that 

= TT^TI / ^(^' OMx)r{k)e~'^''dxdk + 0{h). (68) 
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Hence 



F{x, 0)4'{x)ip*{x)dx. 



The latter equality is obtained by computing the integral over k. This integral 
over k is the inverse Fourier transform of the function 'ip*{k), and it yields 
the function iIj*{x). Lemma 3 is proved. 

The distribution W'lp, as Wigner's distribution, is not nonnegative, and 
the distribution p^, given by expression fl29l) for the wave function ip, is 
nonnegative. 

Denote by Fp^ the average of the function F{x,p) with respect to the 
distribution p. That is, 



^p*=/ F{x,p)p^{x,p)dxdp 

JR?" 



= jF{x,p)i,{y)r{y')x{x - y)x{x - yY^^'~'^^"'dydy'dxdp. (69) 

Lemma 4. IfF{x,p) is a smooth function which, together with its deriva- 
tives, grows at infinity no faster than a polynomial, and ip{x) is an arbitrary 
smooth complex valued function rapidly decreasing at infinity, then the fol- 
lowing equality holds: 

Fw^ = Fp^+0{h) = Fp^+onil), (70) 



where on{l) is an infinitely small quantity with respect to h. 

Proof of Lemma 4. Consider Fp^ given by expression fl69|) . Let us 
represent the function il!*{y'), using composition of the direct and inverse 
Fourier transform f[T^ . in the following form: 

r{y')=J'-'m*iy"),y"^k],k^yl 

i. e., in the following form: 
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Let us substitute this expression into expression (1^ . After simple transfor- 
mations under integral we obtain 

^p^ = /o-^A^r, F{x,p)i){y)i)\y")x{x - y)x{x - y') 



(27r;i)2" Jr6 



xe 



i{y"k- 



-yp+y'p-y'f^)/^dy'dydy"dkdxdp. (71) 



In this expression, let us integrate over y' using formula ( 15^ from Lemma 1. 
We obtain 

F{x,p)i){y)i)\y")x{x - y)x{p - k) 



X e''^y"''~yP+^P"-'''^/^dydy"dkdxdp. (72) 

In this expression, let us make change of variables, introducing new variables 
C, = p — k and r] = x — y. Then, p = k + ^,x = y + ri, and 

xe'^y"^~y^+''^^'^d7]didydy"dk. (73) 

Assuming rj and ^ to be small quantities, let us decompose the function 
Fill + ri,k + ^) by the Taylor formula at the point {y, k) up to the terms of 
second order. We obtain the following expression, in which the values of the 
function F and its derivatives are taken at the point [y, k): 

n 

F{y + rj,k + = Fiy, k) + KA) 

1=1 

1 " 

+^ T.iK,^,v^vJ + 2F:;^,t>M^ + F;:^^^^^,Q + o{v^e,vO■ (74) 

Let us substitute this decomposition instead of the function F{y+ri, k+^) 
into the latter integral, and let us integrate over the variables rj and ^ using 
the formulas written out in Lemma 1. We obtain 



ny,k)+lt{F^^^.;^^+F^^^^^^^^^ 



xip{y)^lj*{y'y^y"'y^''/^dydy"dk 

= Fw' + hSw' + o{h), (75) 
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where S{x,p) = \^(F'^^^,^{a/2b) + iF'^^^^^+F;^^^^{b/2a)) is the second sum- 

mand under the integral in expression (175|l divided by /i, and the average 
values Fw/ and Syj' of the functions F and S with respect to the distribu- 
tion are given by expression flMl) . 

Since by statement of Lemma 4, the function F{x,p) grows at infinity, 
together with its derivatives, no faster than a polynomial, the same property 
holds for the function S{x,p). Let us apply Lemma 3 to the function S{x,p). 
We obtain that Sw^ is bounded as ^ — >^ 0. This and equahty (175|) imply the 
equality Fp^ = F^t^ + 0{h), which is equivalent to the equality required in 
Lemma 4. 

Thus, all the Lemmas necessary for the proof of Statement 2, are proved. 
Let us now proceed to the proof of Statement 2 itself. 

def 

By Statement 1, ctmax = 2-max<^Q | |y4(^o~-Po^<^o| I; where A is given in our 
case by expression (jl]), the projection operator Pq is given by expression (|27|) . 
and normalized functions ipo are given by expression fl2^ . and the function 
(^o<^o is the probability distribution in the physical region of the phase space 
for the given process. 

If the operator A is represented as a sum A = I]|"|^v4j, then, by the 
property of the norm stating that the norm of the sum of vectors is no 
greater than the sum of norms of these vectors, we have 

2n+l 2n+l 

\\A(^o ~ PoA(^o\\ = \\ {A,^o ~ PoAs'^o)\\ < \\As^o - PqA,(^o\\- (76) 

s=l s=l 

Hence, for the estimate of the quantity | \AipQ — Pq^'^oI I it suffices to estimate 
the quantities HA^c^o — -Po^s<^o||, for s = 1, . . . , 2n + 1. 
In our case the operator A is given by expression (j4]): 

^fdHdif dHdip\ z/ ^dH \ 



and A2j-i 



dH d 



dxj dpj 
and A2n+i 



dH d 

dpj dxj 

i 



for j 



n; 



dH 



(78) 
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Note also that since Pq is a self-adjoint projection operator, the vectors 
PoAs<fo and Agt^o ~ Pq^s'^o are orthogonal. The sum of these vectors equals 
Agifo- Hence the following equality holds: 

\\Asipo - PoA^olP = Ps'^olP - ||Po^<^oir. (79) 

Let us start estimating these quantities, starting with s = 2n + 1. 

1. The case A2„+i = -i/hf{x,p). We estimate | |A2n+i<^o--Po^2n+i<^o| P, 
when the operator A2n+i is the operator of multiplication by the function 
—i/hf{x,p) and 

f{x,p) = H-j:—p,. (80) 

j=i ^Pj 

1.1. An estimate of | |y42n+i<^o| In this case, for the quantity 

||^2n+i<^o|P = (^2n+i<^o; ^2n+i<^o), after substitutiou of i/hf{x,p) instead 
of A2n+i, substitution of expressions (12^ for ipQ, and multiplication of both 
parts of the equahty by h'^, we have 

- Ii2 /" J J x2/ N /" , T / X / s i(^-y)p 



Ji \\A2n+i^o\\ =j7:r^7z dxdp f {x,p) dy ilj{y)xix - y)e b 
X / i'*{y')x{x-y')e'^^^^dy' 

dxdpf {x,p) dy i'{y)x{x - y)e—f^Iint, (81) 



where h^t = J^J^ r{y')x{x - y')e^" dy' . (82) 

To transform the integral lint in this expression, let us represent the 
function ip*{y'), using composition of the direct and inverse Fourier trans- 
form ffTSl). in the form 



r{y')=T-'[T[r{y"),y" -^klk-.y'l 

i. e., in the form 
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Let us substitute this expression into expression (152]) for lint- We obtain 

= {^-nUfni^ l^J^yy ' X{x-y')e n dy'dkdy" 

^ ^ r{y")e'^'^ I x{x - y')^^ dy'dkdy". (83) 



Let us compute the latter integral over y' by formula fl54|) from Lemma 1. 
After substitution we obtain 

km = j^J ny")e'^x{p-k)e^dkdy" 

^ ^ riy")xip-k)e-^'^dkdy". (84) 



(27r/;,)" Jr2 

Let us substitute the obtained expression for Imt into expression ( ISTI) for 
/i^l |y42n+i'^o| After simple transformations we obtain 

i{y" k — yp-\-xp — xk) 

xe B (ifc (iy (iy (ix (ip. (85) 

In the obtained expression, let us make change of variables, introducing new 
variables ^ = p — k and rj = x — y. Then p = k + C,,x = y + ri, and 

1 r .2, 



xe (ifc (iy rfy drj dS,. (86) 

Assuming that t] and ^ are small quantities, let us decompose the function 
P{y + r], k + ^) into the Taylor series at the point {y, k) up to terms of the 
second order. We obtain the following expression, in which the values of the 
function / and its derivatives are taken at the point (|/, fc), 

n n 

fiy +v.k+ =/'+E(2/4^.+2//;6)+ E [(/;/; 

'^ifxjfpi "I" f fpi,Xj)^iVj + ifpjfxi "I" f fxi,pj)Vi^j 

Hfpjp. + ff;.jm + o{ri\ e, VO- (87) 
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Let us substitute this expression, instead of function P{y + Ti,k + $,), into 
the latter integral and let us perform integration over the variables rj and ^ 
using the formulas written out in Lemma 1. We obtain 



n ■ II ,1 bh 

~^i.f Xj f pj ~^ f fpj,Xj^^^ ~^ i.fpjfpj ~^ f '^Pj'Pj^ 2,(1 



+ oih) 



x^{y)ij*{y")e ^'"n'"'^ dk dy" dy. (88) 

1.2. An estimate of | |Po^2n+i'^o| Let us now estimate the expression 
subtracted in ([79]), i. e., | |-Po^2n+i'^o| P = {PoMn+iW^ Mn+i^o) , where the 
operator A2n+i is the multiplication by the smooth function —i/hf{x,p). 
Expanding this expression with the scalar product and substituting into it the 
expression for the operator A2n+i, expression for (^q and expression (1271) 
for the operator Pq, represented, in the notations of Lemma 1 in the form 

and multiplying both parts of the equality by /i^, we obtain 

^'||Po^2n+i^oir = (2^^)3l/22n/2 ^„ ^{v)^^^' - y)e'^^ f {x' , p') 

/ X — X'\ ^ fp' — P\ i{x^x'){p+p') 

xip*{y)x{x — y)e » dy dy dx dp dx dp 

= -, — TT — ^ I i'{v)'xW — v)^^ '''^^^ fix',p') 

(27r;i)"2"/2 ifisn ^^^^^^ 

f X — X'\ ^ f P' — P\ i(x-x')(p+p') 

xf(x,p) lint dy dx' dp' dx dp, (89) 



where I,nt = ^^^\y,2 i^*iy')x{x - y')e'"n'''''dy'. 



The integral lint is the same as in fl82|) in the computation of |y42„+iV5o| P- 
Let us substitute into expression fl5^ the representation of the integral lint 
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in the form After simple transformations we obtain, 

/^^||P0^n+1^0|P= (2vr;.pn2n/2 f (x' , p') f (x, pMyW (v") 



^X{x' - - k) 



V2 V2 

i{x — x'){p-\-p') . i{y" k — yp' -\-x' p' — xk) 

xe 2fi + R c^/c c^i/ dy dx dp dx dp. (90) 

In the obtained integral, let us make a change of variables, introducing 
the new variables t] = x' — y, ^ = P — Vj = x — x\ ^ = p' — p. 
Then, x' = y + rj, p = k + ^, x = y + t] + r]', p' = k + ^ + 
After substitution of x', x,p,p', expressed through the new variables, and 
after simple transformations, we obtain 

x^(y)^*(y")x(r7)x (V/v^) X (r/v^) X(0 
xe^+ '^'l'^^'^' +^dfe dy"dy d^ d^ drt'd^'. (91) 

Since the functions and yield the densities of normal distributions 
with small dispersions, let us decompose the function f{y + ri, k + ^ + ^')/(l/ + 
T) + rj',k + ^) in the previous expression into the Taylor series at the point 
{y, k) up to terms of the second order, assuming the quantities 77, rj' , ^, ^' to 
be small. We have, 

f{y + y],k + i + Of{y + v + v\k + o 

n 

= fiy, k) + + + O + ff'^'^J + + ffpM 

+\ E /[/;^,.,^.^.+/;:^,p,^.fe+e;)+/^,.,(e.+e:)^,+/^,,,(6+0(e,+e;)] 

1 " 

+^ E /[i;:^,.,(^^+^D(^.+^;)+i;:^,p,(^^+^Do+/p:,.,e.(^,+r];)+/^,p,6o] 

+0 [r,\ (r^ + V)', (e + e')', ^(^ + r), (r? + V)e) • (92) 
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After substitution of the decomposition of the function f{y+ri, k+$,+^')f{y+ 
rj + r]',k + C,) in the form fl92l) into expression (lUT]) and computing integrals 
over rj, rj', ^, ^' using the integrals of Lemma 2, we obtain 



1 r " 

+ (/;./;, + fJlJ^f^ + ffl,M/{2a)] + o{h)] 
x^{y)^*{y")e^^^'^dk dy"dy. (93) 

Thus, we have estimated the expression |A2n+i<^o| P by formula (IHH!) 
and expression |Po^2n+i<^o| P by formula (193!) . Let us substitute these 
formulas into expression (17^ . let us reduce similar terms, and let us divide 
both parts of the equality by 'h?. We obtain 

||^2n+l'^0 " -Po^2n.+l<^o| P = ||^2n+l'^o|P ~ I |-Po^2n+l'^0 I P 

— — i(y" k — yk) 

xtP{y)tP*{y")e ^ dk dy" dy. (94) 

The last row of expression flM|) is the density of the distribution W'^ 
given by expression fl63|) . Let us apply Lemma 4 to the right hand side of 
the equality flM|) . We obtain 



|^2n+lV'0"~-Po^2n+lV^o| P — T 



(95) 



where p^{,k) is the nonnegative density of distribution given by expres- 
sion (!29l) for the function 

To estimate the required expression max^^ ||^2n+i<^o ~ -Po^2n+i'^o| 1; in- 
troduce the constant C2„+i by the following equality: 



where the maximum is taken over the physical domain U of values of the 
coordinates and momenta for the given process, which contains the support 
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of the density function of the probabihty distribution p^. Then, equahty (^31) 
imphes that 

= {cL^i + on{l))l. (97) 

The latter equahty is obtained from the normahzation condition for the den- 
sity p^p, i. e., from the equahty Jfi„ p^dydk = 1. 

Hence, taking the square root from both parts of the inequahty, we finaUy 
obtain 

\\A2n+l(Po - Po^2n+l'^o|| < (C2n+1 + 0^(1)) (98) 

This finishes examining Case 1. 

2. Case A2j-i = Let us estimate the quantity ||A2j_i(^o ~ 

Po^2j-i<^o|P = ||^2j-i<^o|P — I |-Po^2j-i<^o| P by the same scheme as in the 
previous case: we separately estimate ||A2j_i<^o|P and | |Po^2j-i<^o| 

2.1. An estimate of | |742j_i<^o| Let us substitute into | |742j_i<^o| P = 
{A2j-i'^o] A2j-i^o) expression ([24]) for i^q, and put A2j-i = We have 

dH d f 1 



dx, dp, \{2.hYI- Jn^ - y')e''~^^''dy'j dxdp 

dH ( 1 /■ 7/ N / ^iixi — Vi) i(^'y)p , \ 



i?2" dxj \ (27r/i)"/2 Jj^n fi 

dH d 

X — — ( J,„i) dxdp, (99) 
dxj dpj 

where lint is given by expression fl82l) . Let us substitute here, instead of 
lint, its expression in the form flMl) . After simple transformations and after 
substitution of the derivative of the function x iii the form ([ST]), we obtain 

II . ii9 1 f dH , , .i(xj — yi) i{^-y)p \ 
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dH d _ i(y"-x)k 

>^4'*iy )t: — t; — ixiP " k)) e « dk dy dy dx dp 
OXj opj 



ia 1 f ( dH \ ^ 



hfi^ (27r/i)3"/2 y^jsn \dx 



i'{yW{y"){xj - yj)x{x - y) 



i{y" k — yp-\-xp — xk) 

X (pj — kj)x{p — k)e B dk dy dy dx dp. 

In the obtained expression, let us make the change of variables r] = x — y 
and ^ = p — k. Then, x = y + T],p = k + C,, and 

x(0)x(Oe « ^ ^ dk dy"dy drj d^, 

where in the function (f^) , instead of variables x and p, we have substituted 
y + rj and k + ^, respectively. 

Assuming r] and ^ to be small quantities, let us decompose the function 

(^) (1/ + ^7' ^ + i^to Taylor series up to the zero order. We have 

Let us substitute this expression into the previous one, and integrate it over 
rj and ^, using the integral fl6Tl) from Lemma 1. Finally we obtain 

xe^^^^rffc (100) 

2.2. An estimate of | |Po^2i-i'^o| By construction, this expression 
satisfies the inequalities 

< ||Po^2i-l^o||^ < p2j_l<^oir 

Hence and from relations fl79l) and fllOOl) we obtain 

\\A2j-lipo-PoA2j-l'fo\\'^ = \\A2j-l'fo\\'^-\\PoA2j-l^o\\'^<\\A2j-l^o\\'^ 

xe-^^^dk dy"dy. (101) 
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The last row of expression fllOll) is the density of the distribution given 
by expression fl63l) . Let us apply Lemma 4 to the right hand side of equal- 
ity ffTOTD . We obtain 



I \A2j-l(po- PoA2j-l(po\\'^ = — 



(102) 



where p^{,k) is the nonnegative density of distribution given by expres- 
sion fl2^ for the function ipQ. 

Introduce the constant 6*2^-1 by the following equality: 



.2 def a fdH\^ 



C' =^ max — , (103) 

where the maximum is taken over the physical region f/ of values of coordi- 
nates and momenta for the given process. 

Taking into account this notation, inequality (11021) imphes that 

\\A2j-\(pQ - Po^2i-i'^oir < ^ Cl^-\P^ dydk + 0^(1)) 

The latter equality is obtained, as in the first case, from the normalization 
condition for the density of the probability distribution p^. 

Hence, taking the square root from both parts of the latter inequality, we 
finally obtain 

\\A2j-m - Po^2j-i<^o|| < (C2J-1 + on{l)) (104) 



This finishes examining Case 2. 

3. Case A2j = —^-£r-- Let us estimate the quantity \ \A2j(po—PoA2jifo\\- 

form (1211), i. e., to the function 



If one applies the operator A2j = — f^^r" function ipo of the 



R" 
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then, taking into account formula fl57l) for the derivative of the function Xi 
we obtain the following equalities: 

_ dHd^, OH 1 f 



b 

at. 
b dH 



ah dpj (27r;i)"/2 

R 



ipj dH 



h dp. J (27r;i)"/2^^ 
ib dH d0Q —ipj dH _ 



a dpj dpj h dpj 
That is, A2j(po = ^2j<^o + ^^^o, where 



def ibdHdcpo ^, _ def ipjdH ^ 
a dpj dpj ' ^■^'^^ h dpj ' 



Hence, using the property of the norm that the norm of a sum of vectors is 
no greater than the sum of norms of these vectors, we obtain the inequality 

\\A2j(po - Po^2i<^o| I < I |^2j<^0 - ^0^2j<^0| I + WKj'pQ " ^0^2i<^o| |- (105) 

Note that the operator coincides with the operator ^2^-1 (see the 
previous Case) if one replaces the function |^ in it to Therefore, 
using formulas fll04p and (11031) . we obtain 



|^^^.<^o - Po4,<^o|| < + 0,(1)) /v^, (106) 

where [CL)^ =^ max 

^ 2a{x,v)&\dpj) 



On the other hand, the operator coincides with the operator Am+i 
see Case 1), in which f{x,p) = Pj§^- Hence, using formulas fl98|) and fl96|l . 



we obtain 
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max 



n 

E 



d / dl£\ya_ /_d_/ dH\Y b 



Put C2j =^ G^j + C'ij. Then using inequalities ([IUSD, (HUn}, and (fTOTj) . 
and using the introduced notation C2j, we finally obtain 

(C'+on{l)) (C''+on{l)) 1 

This finishes examining Case 3. 

Now we are ready to finish the proof of Statement 2. 

def 

By Statement 1, amax = 2 ■ maXi^^ H^^o ~ -Po^'^oll- Hence, applying to 
inequahty ([75]) relations dM]), ffTIHD . and ffTUg]) . we obtain 

2n+l 

a^ax = 2 ■ m_ax | |A(^o - ^o^'^ol I < 2 ■ m_ax V \\A,ipQ - Pq^^'^oI I 

s=\ 

2n+l 1 1 

= 2j:iCs + 0,(1)) ^ = iC + 0.(1)) ^, 



s=l 



where C =^2Es=t^Cs. 

Q. E. D. Statement 2 is proved. 



38 



Appendix 2 



Proof of Theorem 5 

For proof of Theorem 5, consider formula fl33l) from Theorem 4 for the case 

2 

when H{x^p) = ^ + V{x). We have 

xx{x - y)x{x - y')e^^y~y''^^i){y\ t)dy'dxdp. (109) 

Thus, in this case Hijj is represented as the sum of three integrals H'tp = 
h + h + h, where 

= 72^ I ^xix-y)xix-y')er^^'-''^P^/j{y',t)dy'dxdp, (110) 

ib 



{2'KhY 



I —Y.Pk{xk-y'k) 

J am 

,Q„ fe=i 



yix{x-y)x{x-y')e^'^y-y>i^{y\t)dy'dxdp, (111) 

xx(x - y)x(x - y')eT^^y-y'^^,l,{y\ t)dy'dxdp. (112) 

Note that the expression /^jn x(x — y)x{x — y')dx from the first integral 
can be transformed to the following form: 

" , [ b r _H--y)^ _K--y'? 

X\x — y)x\x — y )dx = I — - 1 J e 2an e ^an dx 

Jin ^ ^ ij" 

b \"^^/" b{a:-(y+y')/2)'^ b(y-y')2 ( V + y'\ b{y-y')'^ 

' ^ e — ^^d[x-- — -]=e 2^.(113) 



Hence the integral l\ is transformed to the form 

= (2^ / |^e-^ei(^-^')^^(y',t)rfy'rfp. (114) 
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As known from the formulas for Fourier transform, for any smooth func- 
tion f{y') the following equality holds: 

^(^) = T^Tt^ I e"^^'-''^^f{y')dy'dp. (115) 

^^^^^ R2n 

Below we shall often use this equality. 
In particular, one has the equality 

'^(y^t) = 7^ I e-'^e^n(y-y')p^(y',t)dy'dp. (116) 

Let us differentiate both parts of the latter equality with respect to yk- 
We obtain, taking into account relation flllSp . 

- ' f^'P" ^^y'-y''\-'^e^^y-y'^^Hy',t)dy'dp 



dyk {2'KhY V h 2ah 

ipk Hv-v'f i 



-e ■4"'* e 



Uy-y')p^{y\t)dy'dp. (117) 



If we differentiate both parts of the obtained equality with respect to y^ once 
more, then we obtain the equality 

d^^iy.t) _ 1 fipkfipk Kvk-y'k 



dyl {2'KhY J^^ n \ h 2ah 

2m 

X( 



ce-'-^^e^^y-y'^Pilj{y',t)dy'dp = -^Ii - (118) 

/ 



where 



^'=(2^ / ^e-'^e^niy-y')p^p(y',t)dy'dp, (119) 

and 



To compute the integrals J^, consider the following equality which is a 
particular case of equality flllSp : 
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Let us differentiate both parts of this equahty with respect to Uk- We obtain 



(2TThy^ J \2ah 2ah \ k 2ah 

X e iah e 



Uy-y')p^{y'^t)dy'dp. (122) 



Hence, taking into account notation (11201) for the integral and relation (11151) . 
we obtain the equalities 

Substituting this equality into relation (11181) . we obtain 
Let us express from this equality: 

_ d^i^ ^ bh 

^ 2m dyl Ama ' 

Taking the sum of the obtained equality over all k from 1 to n, we obtain an 
expression for the integral Ii. 

Let us pass to computing the integral I2 given by expression (II 111) . 
Consider equality (I113P of the form 

^ y)X{^ ~ y )dx = e 4aR ^ where X[x^y)=\ — t\ ^ ■ 
J \ ann / 

Let us differentiate both parts of this equality with respect to y'j^. We obtain 
7 an 



2ah 
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or, after omitting common factors, 



>fe - yk)x{x - y)x{x - y')dx = ^^ e '"4-^ ' . (125) 



Substituting this equality into expression fillip for the integral /2, we 
obtain 



(27r^) 



Comparing the latter expression with expression ( 1120p for the integrals 
we obtain I2 = — ^^/?^Z]fc=i -^4 • Substituting here the computed expres- 
sions (11231) for the integrals J4, we finally obtain 

nbh , , , , 

h = (126) 

Now consider the integral I3 given by expression (11120 . In accordance 
with formula (11150 . this integral can be transformed to the form 

^3 = ^{y,t) J iVix) - J2 'o^^^k ~ yk)]x^{x - y)dx, 

where 'X^{x — y) = {h/ {aTih)Y^'^exp{—h{x — / {ah)) is the density of the 
normal distribution, or, after the change of variables x'f. = Xk — yk, to the 
form 

h = i^{y,t) J {V{y + x') - ± 



k=l 



Whereas the previous integrals have been computed exactly, let us com- 
pute this integral approximately, assuming the dispersion ah /2b of the normal 
distribution to be a small quantity, decomposing the function V{y+x') 

into the Taylor series at the point y with respect to x' up to the second order, 
and decomposing dV{y + x')/dyk up to the first order. We have 

l„ \k=i k=ik'=i^yk^yk' J 
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Hence, since for the density of the normal distribution x^i^') the following 
relations hold: J^n x^{.x')dx' = 1, x'f^x^{x')dx' = 0, x'i,x'i.,x'^{x')dx' = 
for k 7^ k', and /r„ x'f,x'f,x'^{x')dx' = ah/{2b), we finally obtain 

h - ^(y, - i^{y. t)^j: (127) 

Thus, since Hip = Ii + I2 + I3, where the integrals /i,/2,/3 are given 
respectively by expressions flllOp . fillip . (11121) . then, substituting here their 
computed values as expressions (I124p , (11260 , (I127P and reducing similar sum- 
mands, we obtain the expression for H required in Theorem 5: 

^2 " ^2 . ah^d^V 3nbh 
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Appendix 3 



An estimate of parameters of the model 

1. An estimate of the parameter a/h of the 
model 

The difference between the operator H given by expression (15^ of Theorem 5 
and the standard Hamilton operator in the Schrodinger equation, can cause 
difference between the spectra of the energy operators (in particular, for the 
potential energy function of hydrogen atom) for the model considered above 
and for the standard model of quantum mechanics. The difference between 
these operators is in the third summand with the factor ah/ Ah. If one assumes 
that behavior of particles is described by the Schrodinger equation with the 
operator H of the form fl5B]) more exactly than with he Hamilton opera- 
tor, then experiments should show the non-exactness of spectra computed 
by means of the Hamilton operator, i. e., non-exactness of non-relativistic 
quantum mechanics. 

The discrepancy between theoretical and experimental data in non-relat- 
ivistic quantum mechanics is well known. It has been discovered in 40s in 
|10j . and has been called the Lamb shift of levels of hydrogen atom. Later 
on, this effect has been explained in quantum electrodynamics by interaction 
of the electron with fluctuating electromagnetic field (see, for example, |16], 
p. 593, where one can find references to original works). 

If one assumes that the data in the Lamb experiment are related with 
the perturbing summand in the operator f/", then these data allow one to 
estimate the quantity ah/Ab. 

In this section the computations of the estimate of the quantity ah /4b 
have been carried out following the computations given in [18] for substanti- 
ation of the size of the Lamb shift in the spectrum of hydrogen atom. 

Let 



e2 



be the potential function of hydrogen atom, and H = E + V he the operator 
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from Theorem 5, where 

Consider the operator H a.s a. perturbation of the Hamilton operator 
H = E + V of hydrogen atom. Let us estimate eigenvalues of the operator 
H. 

The standard perturbation theory implies that at the first approximation, 
the correction 6En to the eigenvalue En of the Hamilton operator has the form 

SEn= I pn{y){V{y)-V{y))dy, 

where Pn{y) = \'^n{y)\'^ and ipniu) is the eigenfunction of the Hamilton oper- 
ator with the eigenvalue En- 

Substituting the expression for V{y) into the expression for 6 En, we ob- 
tain 

ah f , ^ d'Vjy) ^ 

Since the integral of the Laplace operator of the function V{y) = —e^/r 
equals 47re^(5o(?/), where 47re^(5o(?/) is the delta function at zero, this implies 
that 

^ ^ ahne^ , , 

5 En = 7—Pn{0). 



For hydrogen atom it is known (see, for example, ^20j, p. 342) that 

nn-^ \ h J 

Hence 

afrne^ 1 ( me^\^ am^e^ am^a^c^ 

oEn = ■ 



b 7m^ \ I b n^h^ b rfih 



where a = e^/ihc) = 1/137 is the fine structure constant, c is the velocity of 
light. Hence, a/6 is expressed as follows: 

lOEn 
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In the Lamb-Retherford experiments for hydrogen atom it has been estab- 
hshed that SE2 = 1058MHz = 1058 ■ lO^h erg, where h = 27ch. Comparing 
this value with the obtained value of 6E2, we obtain by simple calculations 
the estimate of the quantity a/b = 3,41 ■ lO'^sec/g. Hence the standard de- 
viation for the density of normal probability distribution Xii^')^ with which 
the smoothing of the potential V is made, equals ^ah/2b = 4, 24 ■ lO^^^cm. 
This quantity is much less than the radius of hydrogen atom. 

2. An estimate of the diffusion coefficients and 
of the time of the transformation process 

In this section we shall assume that the diffusion coefficients a and b are 
defined by the standard heat action of the surrounding medium on the moving 
electron. One can assume that the Brownian particle (the electron) is acted 
on by a fiuctuating force from the surrounding medium, and also by stochastic 
resistance proportional to the velocity of the particle. For modelling of the 
motion in this situation, one usually uses (see, for example, [T7], p. 196) the 
Langevene equation 

p = --fp + F{t), 

where 7 is the friction coefficient for the unit mass, and F{t) represents the 
fluctuating force, which is assumed to be independent of the velocity, with 
the mean value equal to zero. 

This equation is equivalent (see [T7|, p. 212) to the Fokker-Planck equa- 
tion for the density of probability distribution f{p,q,t) in the phase space, 
in its standard form: 

dt m dxi dpi dpi ' 

where m is the mass of the particle (in the case of electron, m = 9.10939 • 
IQ-^^kg), k is the Boltzmann constant {k = 1.38066 ■ 10~^^J/K), T is the 
temperature of the medium. 

The solution of these equations for time intervals t much greater than 7"^, 
is well known (see |2Tj, p. 215) to yield the diffusion process with respect to 
coordinates x with the diffusion coefficient kT/{m'~f). This relation is also 
known as the Einstein relation for the diffusion coefficient ([I7j, p. 198). 
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Thus, under these assumptions one can suppose that the diffusion coeffi- 
cient with respect to coordinates = kT/{mj), and the diffusion coefficient 
with respect to momenta (as in the Fokker-Pianck equation) = 'jkTm. 

Hence, a/b = {jm)^^ and ab = kT. 

By the estimate obtained in the previous section, we have 

y = — = 3.41 • lO'^sec/kg, 

b jm 



and the quantity 
1 

7 



3.41 -W -m 3.41 • 10^ • 9.10939 • IQ-^i 



= 3.22 • W^^sec-^. 



On the other hand, the time of the transformation process to the process 
described by the Schrodinger equation, is expressed, by Theorem 4, by the 
quantity h/{ab) = h/{kT). This time for T = 1°K equals 7.638 ■ lO^^^sec. 

For the same temperature, the diffusion coefficients are estimated in the 
following way: 

= kT/{m-f) = 4.708 • IQ-^^m^ • sec"^; 
6^ = -ikTm = 4.049 • lQ~^^{kg ■ m ■ sec-'^f ■ sec~^. 
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